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Abstract. We use a Press-Schechter-like calculation to study how the 
abundance of voids changes in models with non-Gaussian initial conditions. While 
a positive skewness increases the cluster abundance, a negative skewness does the 
same for the void abundance. We determine the dependence of the void abundance 
on the non-Gaussianity parameter f a \ for the local-model bispectrum — which 
approximates the bispectrum in some multi-field inflation models — and for the 
equilateral bispectrum, which approximates the bispectrum in e.g. string-inspired 
DBI models of inflation. We show that the void abundance in large-scale-structure 
surveys currently being considered should probe values as small as f n \ < 10 and 
/_? < 30, over distance scales ~ 10 Mpc. 
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1. Introduction 

The paradigm of cosmological structure formation from a spectrum of primordial 
perturbations like those predicted by inflation has now been fairly well established by 
cosmic microwave background (CMB) experiments [JJ. We are thus now motivated 
to test more precisely the predictions of inflation and to look for possible deviations. 
One of several such possibilities is measurement of departures from Gaussianity of the 
initial perturbations (see, e.g., Ref. [2] and references therein). The simplest slow-roll 
single-field models of inflation predict that primordial perturbations should be very 
closely Gaussian [3] , but with predictably small departures from Gaussianity [H [5] . 
Multi-field [6 J models, such as the curvaton model [7J, and string- inspired DB1 8 
inflationary models can produce larger deviations from non-Gaussianity. 

Departures from primordial Gaussianity can be sought in the CMB [9] [JO] [JT] , 
large-scale structure (LSS) p~2], and the abundances and properties of the most massive 
gravitationally-bound objects in the Universe today or at high redshift [H] HH [15] . 
The CMB provides a more powerful and clean probe of primordial non-Gaussianity 
than direct measurement of the bispectrum in low-redshift LSS surveys in models with 
scale-invariant non-Gaussianity [10], although biasing may amplify the effects of non- 
Gaussianity on LSS to the level where they may be comparable in detectability to the 
CMB [TBI H3 HH1 US] • Measurements of the cluster abundance may do better than the 
CMB and LSS if the non-Gaussianity is not scale-invariant [50], as may occur in DBI 
models. What is clear, however, is that the thorny systematic effects that enter in all 
of these approaches will require that a variety of complementary avenues be taken to 
establish a robust detection of non-Gaussianity. 

Voids have been considered as probes of cosmology, but no systematic study has 
been carried out for voids as probes of primordial non-gaussianity |21j . In this paper, 
we consider the abundance of voids as a test of the distribution of the primordial 
perturbations. Galaxy clusters form at the highest overdensities of the primordial 
density field and thus probe the high-density tail of the primordial density distribution 
function. Similarly, voids form in low-density regions and should thus probe the low- 
density tail of the distribution function. If there is a large negative skewness, the 
void-size distribution function will be increased at the largest void sizes and decreased 
at smaller void sizes, opposite to the effect on the cluster mass function. 

In Section[2] we develop a Press-Schechter (PS) estimate of the void abundance for 
Gaussian initial conditions. This PS-like calculation is far from state of the art [22] [23] 
for Gaussian initial conditions. However, it is easily generalized to non-Gaussian 
initial conditions and should be sufficiently reliable to estimate the fractional effects 
of non-Gaussianity on the void abundance (as it describes well the halo abundance 
[24] , In Section[3l we discuss the relation between the skewness and the non-Gaussian 
parameter f n \ for the local model [4], which approximates the non-Gaussianity in 
multi-field models, and the equilateral model [5j, which approximates that in e.g. 
string-inspired DBI models. In Section [4] we provide results of the void-abundance 
calculation, and we estimate the smallest f n \, for the local model and for the equilateral 
model (with and without scale dependence) that should be detectable in several 
surveys currently under study. In Section [5j we make some concluding remarks and 
outline further steps that must be taken before the void abundance can be used to 
probe non-Gaussianity. 
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2. The Press-Schechter Abundance 



We begin by developing a PS-like calculation of the void abundance, but we first 
review the standard Press-Schechter calculation of the halo mass function. 

Formation of a bound halo requires a linear-theory density fluctuation 5r > S c ~ 
1.69, smoothed on scale R. The smoothing radius R defines the halo mass M via 
M = (47r / '3) pbR 3 , where p& is the mean nonrelativistic- matter density. 

The differential abundance of dark-matter halos as a function of mass and rcdshift 



is 



dn _ £b_ 
dM ~ 1 M 



dP(> S c \z,M) 



dM 



(1) 



where P(> 5 c \z, M) denotes the probability that 5r lie above the threshold for collapse 
S c . For Gaussian initial conditions, 



P(> S c \z,R) 



1 or 

/2tt 5c 



exp 



24 



(2) 



where <tr denotes the rms mass fluctuation on a scale R 1 and there is an implicit 
redshift dependence in <tr. We then introduce the Press-Schechter swindle — i.e., 
that every mass element in an underdense region gets absorbed into the nearest 
overdensity — by introducing a factor of 2 in the abundance. The differential abundance 
of halos is then 



dn 



-AM 



d\a. a M 



dlnM 
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The mass function is then normalized so that 
mass element in the Universe is housed somewhere 



e - & c/^MdM. 
M(dn/dM)dM 
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p b ; i.e., every 



2.1. Void abundance: Gaussian initial conditions 

A similar calculation can be used to estimate the void abdundance. In the PS 
description, every mass element is in a gravitationally bound structure of some mass, 
and the densities of these objects are all > 200 times the mean density. This means 
that > 99.5% of the volume of the Universe is empty. A void distribution function 
can be derived in a way analogous to the PS mass function by realizing that negative 
density fluctuations grow into voids (as opposed to positive-density fluctuations, which 
grow into bound objects), and that there is a critical underdensity S v for producing a 
void, to replace S c , the critical overdensity for producing bound objects. In principle, 
S v can be calculated from theory, but for now, we will treat it as a phenomenological 
parameter. 

The radius R of a spherical volume in which a mass M has been cleared out is 
R = (3M/47rpfc) 1 / 3 . We thus derive (dn/dR), the differential abundance of voids of 
diameter i?, 



dn 
dR 



9 
2^2 



7T 1 \6y\ 
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din aM 
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dlnM 

We here introduce again the Press-Schechter swindle, to take into account the fact that 
each underdense region expands by a factor of two so that the entire (actually, 99.5% 
of the) volume of the Universe is occupied by a voids of radii R with the distribution 
dn/dR; i.e., the size distribution is normalized so that 
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Figure 1. Skewness S^ij for a non-Gaussianity parameter / n i = 1 as function 
of radius R. The dashed line is the local non-Gaussian model; the solid line is 
the equilateral model; and the dotted line is the DBI-type equilateral model for a 
scale-dependence parameter re = —0.3. See Ref. |20| for details. We also indicate 
the scales probed by halo and void abundances. 



2.2. Void abundance with 



-Gaussian initial conditions 



The modification of Eq. (j4| when there is a small primordial skewness S3 m (which 
may most generally depend on the mass scale M) then follows from the analogous 
modification for the halo abundance in Refs. [TH [15l [20] . The only subtlety is that S v 
is now a negative quantity (since clusters come from overdensities, while voids come 
from underdensities). In fact P <s = 1 - P>s thus \dP <s JdM\ = \dP >s JdM\. We 
can still use \5 V \ provided that Eq. fH is replaced by (see [20] ) 
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3. The skewness 



According to Eq. ([H]), the void abundance depends on the non-Gaussianity only 
through the skewness 5*3^. The skewness then depends on the detailed form of the 
bispectrum. There are a variety of bispectra considered in the literature. The standard 
"working-horse" model for non-Gaussianity is perhaps the local model, which features 
a bispectrum that arises from multi-field inflation models. There is then the equilateral 
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model, which approximates the bispectrum in DBI models. DBI models can also allow 
for a scale-dependent non-Gaussianity, with power-law index k. See Ref. [5J for 
definitions and further details. 

In Fig. [TJ we plot the skewness S^ij for a non-Gaussianity parameter f n \ = 1, 
as a function of R for the local model, the equilateral model, and the equilateral 
model with a scale-dependence parameter k = —0.3. The skewness in the equilateral 
model is about a factor 3 smaller than that for the local model. We thus infer that 
constraints to / n i for the local model that derive simply from S3. a will be a factor 
3 more stringent than for the equilateral model. An equilateral model with scale 
dependence characterized by n = —0.3 will have constraints similar to those for the 
local model. 

We should also note that \6 V \ is expected to be smaller than 6 C by a factor 2 — 3. 
If voids and halos probe the same scales, then voids should yield constraints on f n \ 
three times worse. However, voids may probe slightly larger scales than halos, and as 
£3 in these models increases with scale, this compensates for |<J„| being smaller than 

Sc- 

4. Results 

To illustrate the promise held by the void abundance for constraining f n \, we consider 
several large-volume LSS surveys at high redshift that are now being considered. 
The curves in Fig. [2] show the ratio of the void abundance with non-Gaussianity 
to the abundance with Gaussian initial conditions, as a function of void size i?, for 
representative values / n i = +30 and f n \ — —30. The curves are evaluated for a 
redshift z = 0.8. The curve with greater void abundance at larger R and smaller void 
abundance at smaller R is for negative f n \ and negative skewness in the density field, 
while that with smaller (larger) void abundance at larger (smaller) R is for positive 
/ni (positive skewness in the density field). The crossover in R between enhancement 
and suppression of the void abundance reflects the onset of the exponential tail of the 
PS distribution. 

The points with error bars estimate the la error bars (coming from Poisson errors 
in the void abundance) anticipated for a survey with a redshift slice Az = 0.3 centered 
at z = 0.8 that covers 30,000 square degrees of the sky. The curves are evaluated for 
a critical underdensity 5 V = —0.7. As indicated by the Figure, signal to noise in the 
determination of f n \ will come primarily from the lower- R end of the distribution, 
where the void abundance is larger. 

Fig.[3]shows the smallest local-model / n i detectable at the la level as a function of 
the critical underdensity S v for several surveys currently under study. The solid upper 
line is a survey with parameters comparable to those proposed for the BOSS SDSS- 
3 survey [27 (which we parametrize as a 10,000 square-degree survey complete to 
redshift z — 0.8); this curve is evaluated using the void-size distribution over the range 
2 h^ 1 Mpc < R < 60 h^ 1 Mpc. The lower two (dashed) lines assume parameters of a 
survey like ADEPT [28] (which we parametrize as a 30,000 square-degree survey over 
the redshift range 1 < z < 2. The upper dashed curve uses only the void distribution 
with R> Sh^ 1 Mpc; the lower dashed curve uses 2 h^ 1 Mpc < R < 60 h^ 1 Mpc. 

Since the void abundance depends on non-Gaussianity only through the 
dependence on S^^r, we can use the curves shown in Fig.Q]to compute constraints for 
the equilateral model, both with and without scale dependence. As expected, since 
the equilateral model (with no scale dependence) predicts a value of 63^ that is about 
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Figure 2. The ratio of the void abundance with non-Gaussianity at the level of 
/ n l | = 30 to the void abundance with Gaussian initial conditions as a function of 
the void size R. The curves are evaluated for a central redshift z = 0.8, and the 
points are la Poisson errors for a survey of width Az = 0.3 that covers 30,000 
square degrees for a fiducial Gaussian case. The curve with a large ratio at large 
R is for / n i < 0, while the curve with the small ratio at large R is for f a \ > 0. 



3 times smaller than in the local model, the smallest the values of accessible with 
the void abundance is larger by roughly a factor of 3 than in the local model. Likewise, 
Fig.[3]shows that f n \ for the equilateral model with n = —0.3 is comparable to that in 
the local model; the constraints to f c , q for k = —0.3 are similar to those in the local 

" nl 

model. This is particularly interesting as the large-scale bias effect due to equilateral 
non-gaussianity is orders of magnitude smaller than the effect for the local case |26) . 
Thus the combination of the two measurements will help discriminate among different 
type of non-gaussian initial conditions. 

As Fig. [3] shows, the results do not depend very sensitively on the precise value of 
S v . The value that we choose for S v depends on the precise definition of a void. Here 
we have used as a canonical value S v = —0.7 for the critical linear-theory underdensity 
for a void. At this value of S v , the physical void is underdense by a factor of 2 (see, 
e.g., Fig. 1 in Ref. [22]), which we believe to be a conservative definition of a void. If 
voids are defined in the survey to be regions that are even emptier, then the relevant 
magnitude \S V \ will be even larger, and according to Fig. [3l the smallest detectable 
| / n i| will be even a bit smaller. 

5. Discussion 

The bottom line of our analysis is that the abundance of voids in LSS surveys that 
are currently being considered can probe values of the non-Gaussianity parameter 
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Figure 3. The smallest local-model / n j detectable at the la level, as a function 
of the critical underdensity S v , for several surveys currently under study. The 
solid upper line is for something like the BOSS SDSS-3 survey using the void- 
size distribution over the range 2 h —1 Mpc < R < 60 h~ 1 Mpc. The lower two 
(dashed) curves are for an ADEPT-like survey. The upper dashed curve uses 
only the void distribution with R > 8/i _1 Mpc; the lower dashed curve uses 
2 h- 1 Mpc > R < 60 h' 1 Mpc. 



in the local model as small as / n i ~ 10 and in the equilateral model (with no scale 
dependence) as small /f? ~ 30. This probe of non-Gaussianity may be competitive 
with those coming from the CMB, LSS, and cluster abundances [TDJ [TSl [2U]. They 
will complement CMB constraints by (a) providing a different avenue with different 
systematic effects, and (b) probing non-Gaussianity primarily over distance scales 2— 
60 Mpc, scales generally smaller than those that will be probed by the CMB. The 
fact that voids constraints on the different types of non gaussianity are comparable 
is particularly interesting. In fact the large-scale bias effect due to equilateral non- 
gaussianity is orders of magnitude smaller than the effect for the local case. Thus the 
combination of the two measurements will help discriminate among different type of 
non-gaussian initial conditions. 

Finally, we note that we have here done no more than estimate the smallest / n i 
detectable with the void abundance. To do so, we have taken a simple Press-Schechter 
approach to estimate the fractional change in the void abundance. While this approach 
should provide reasonable rough estimates to the fractional change in the abundance, it 
will not reliably provide the abundances themselves. Before the void-abundance probe 
of / n i can be implemented, we will therefore require significantly more sophisticated 
calculations, which will probably require numerical simulations, to accurately model 
not only the void growth, but also the systematic effects that will arise in any realistic 
void-identification algorithm. We hope that our results motivates this type of future 
work. 
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